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ABSTRACT 
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^ • Context. So far, red giant oscillations have been studied from radial velocity and/or light curve variations, which reveal frequencies of the 
^N) oscillation modes. To characterise radial and non-radial oscillations, line profile variations are a valuable diagnostic. Here we present for the 

■ first time a line profile analysis of pulsating red giants, taking into account the small line profile variations and the predicted short damping and 

■ re-excitation times. We do so by modelling the time variations in the cross correlation profiles in terms of oscillation theory. 
00 Aims. The performance of existing diagnostics for mode identification is investigated for known oscillating giants which have very small line 

profile variations. We modify these diagnostics, perform simulations, and characterise the radial and non-radial modes detected in the cross 
correlation profiles. 

Methods. Moments and line bisectors are computed and analysed for four giants. The robustness of the discriminant of the moments against 

■ small oscillations with finite lifetimes is investigated. In addition, line profiles are simulated with short damping and re-excitation times and 
Oh their line shapes are compared with the observations. 

^ Results. For three stars, we find evidence for the presence of non-radial pulsation modes, while for ^ Hydrae perhaps only radial modes are 
5_j ' present. Furthermore the line bisectors are not able to distinguish between different pulsation modes and are an insufficient diagnostic to 
' discriminate small line profile variations due to oscillations from exoplanet motion. 

a '. 
> ' 

• i-H Key words, stars: oscillations - stars: individual: e Ophiuchi, r] Serpentis, <f Hydrae, S Eridani - Techniques: spectroscopic - Line: profiles 

■ - - 1. Introduction For several stars on or close to the main sequ ence, solar like 

oscillations were detected for a decade (see e.g.'Kieldsen et al.l 
Techniques to perform very accurate radial velocity observa- 0995), Bouchv & Carrier (2001), Bouchv & Carrier (2003), 
tions are refined during the last decade. Observations with Bedding & KieldsenI (l2003l) ). More recendy, such type of 
an accuracy of a few i ns'' are obt ained regularly (see e.g. oscillations has also been firmly established in several 
[Marcv & Butler (2000), lOueloz et al. (2001b)) and detections red (sub )giant stars . Frandsen et al. (2002), De Ridder et al.j 
of amplitude s of 1 ms" ' are nowadays possible with e.g. (200 6bl). ICarrier et alJ tOOdti . (see also iBarban et alJ ( 12004 1 
HARPS (Peoe et alJl2003l). This refinement not only forced a and Carrier et al.1 ( l2003l) used the CORALIE and ELODIE 
breakthrough in the detection of extra solar planets but also in spectrographs to obtain long term high resolution time series 
the observation of solar like oscillations in distant stars. Solar of the three red giants f Hydrae, e Ophiuchi and 77 Serpentis 
like oscillations are excited by turbulent convection near the and of the subgiant 6 Eridani, respectively. They unravelled a 
surface of cool stars of spectral type F, G, K or M and show large frequency separation in the radial velocity Fourier trans- 
radial velocity variations with amplitudes of typically a few form, with a typical value expected for solar like oscillations 
ems ' to a few ms ', and with periods ranging from a few in the respective type of star (a few juHz for giants and a 
minutes for main sequence stars to about half an hour for sub- few tens fiUz for s ubgiants), according to theoretical predic- 
giants and a couple of hours for giants. tio ns (e.g.LO ziembowski et alJ ( l200ll) ). It was already predicted 

b V 'P ziembowskii 1 1 97 L) that non-radial oscillations are highly 
damped in evolved stars, and that, most likely, any detectable 
oscillations will be radial modes. 
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So far, red (sub)giant oscillations have only been stud- 
ied from radial velocity or light variations. Line profile vari- 
ations are a very valuable diagnostic to detect both ra- 
dial and non-r adial heat driven coherent oscillations (e.g. 
lAerts & De Cad (|2003) and references therein), and to char- 
acterise the wavenumbers (£,m) of such self-excited oscilla- 
tions. It is therefore worthwhile to try and detect them for red 
(sub)giants with confirmed oscillations and, if successful, to 
use them for empirical mode identification. This would pro- 
vide an independent test for the theoretical modelUng of the 
frequency spectrum. 

It is also interesting to compare the line profile diagnostics 
used for stellar oscillation analysis with those usually adopted 
to discriminate oscillations from exoplanet signatures, such 
as the line bisector and its derived quantities. Recently, iGravl 
pointed out that the wide range of bisector shapes he 
found must contain information about the velocity fields in the 
atmospheres of cool stars, but that the extraction of information 
about the velocity variations requires detailed modelling. Here 
we perform such modelling in terms of non-radial oscillation 
theory. We do so by considering different types of line char- 
acteris tics derived from cross correlation profiles. IPall et af] 
already concluded that bisectors are not suitable to anal- 
yse solar like oscillations. We confirm this finding and propose 
much more suitable diagnostics. We nevertheless investigate 
how line bisector quantities behave for confirmed oscillators, in 
order to help future planet hunters in discriminating the cause 
of small line profile variations in their data. 

The main problems in characterising the oscillation modes 
of red (sub)giants are, first, the low amplitudes of the ve- 
locity variations, which results in very small changes in the 
line profile. Secondly, the damping and re-excitation ti mes 
are predicted to be very short. Indeed, IStello et aP ( l2004l) de- 
rived an oscillation mode lifetime in ^ Hydrae of only ap- 
proximately two days. In general, mode lifetimes are diffi- 
cult to compute for different stellar param eters f rom current 
available theory, cf. Houd ek et al versus Stello et al. 

12004 ). Stello et al , (2006) derived the observed mode lifetime 
of ^ Hydrae from extensive simulations and find a large dif- 
ference with theoretical predictions. The known mode iden- 
tification methods from line profile variations all use an infi- 
nite mode lifetime so far Here, we provide fine profile vari- 
ations simulated for stochastically excited modes, and inves- 
tigate how the damping affects the behaviour of the diag- 
nostics in this case. We apply our methodology to the four 
case studies of ^ Hydrae (HD 100407, G7III), e Ophiuchi 
(HD146791, G9.5III), 77 Serpentis (HD168723, KOIII) and 6 
Eridani (HD23249, KOIV). Some basic properties of these stars 
are listed in Table[2 

The paper is organised as follows. In Sect. 2 the data set 
at our disposal and different observational diagnostics are de- 
scribed. In Sect. 3 the diagnostics for mode identification are 
described, while in Sect. 4 the damping and re-excitation ef- 
fects in theoretically generated spectral lines are investigated. 
In Sect. 5 the results obtained from the simulations and ob- 
servations are compared and in Sect. 6 some conclusions are 
drawn. 
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Fig. 1. Pow er spectrum of e Ophiuch i. The lower black one is 
obtained bv IPe Ridder et all ll2006hl) . the grey one in the mid- 
dle is obtained from (v) and the top one is obtained from the 
bisector velocity span. For clarity the latter two are shifted. 



2. Observational diagnostics 

2.1. Spectra 

For all four stars we have spectra at our disposal obtained with 
the fibre fed echelle spectrograph CORALIE mounted on the 
Swiss 1.2 m Euler telescope at La Silla (ESO, Chile). The ob- 
servations for ^ Hydrae were made during one full month (2002 
February 18 - March 18). This is the same dataset as used by 
Frandsen et al. (2002) for the detection of the solar like oscil- 
lations. For e Ophi uchi and n Ser pentis the sol ar like obser- 
vations described by De Ridder et al.. C2006b.) and lCarrier et alJ 
(2006) are obtained from a bi-site campaign, using CORALIE 
and ELODIE (the fibre fed echelle spectrograph mounted on 
the French 1.93 m telescope at the Observatoire the Haute- 
Provence) during the summer of 2003. The observations of 
6 Eridani are obtained with CORALIE during a twelve day 
campaign in November 2001. For the line profile analysis, de- 
scribed in this paper, only the data obtained with CORALIE are 
used. These spectra range in wavelength from 387.5 nm to 682 
nm and the observation times were adjusted to reach a signal to 
noise ratio of at least 100 at 550 nm (60-120 for S Eridani) with- 
out averaging out a too large fraction of the pulsation phase. 
More details about the observation strategy is available in the 
publications describing the first detections of the solar like os- 
cillations of the respective stars. 



2.2. Cross correlation profiles 

The line profile variations o f the four pulsatin g red (sub)giants 
are analy sed wi th moments (l Aerts etalJljgl) and with the line 
bisector ( lGravl l2005). The moments are often used to analyse 
stellar oscillations excited by a heat mechanism, while the line 
bisector is usually adopted to discover planet signatures. Both 
methods can be applied to one single spectral line, but also to a 
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Table 1. Basic stellar parameters of the four stars: Effective temperature (Tes) in Kelvin, rotational velocity (i/sin/) in 
km s parallax (n) in mas, distance (d) in pc, the apparent magnitude (;«,,) and absolute magnitude (My) in the V band. 



parameter 



6 Ophiuchi" 



Teff [K] 

usini [km s" 
n [mas] 
d[pc] 
m„ [mag] 
My [mag] 



4887 : 
3.4 4 
30.34 : 

33.0: 

3.24 4 
0.65 4 



: 100 

0.5 

:0.79 
:0.9 

0.02 
0.06 



Serpentis 
W 



4855 
2.6 H 

52.81 
18.9 

3.23 = 

1.85: 



0.8"* 
: 0.75'' 
tO.3 
0.02" 

:0.05 



^ Hydrae 
5010 + 15* 
2.4'' 
25.23 + 0.83" 
39.6 ± 1.5 
3.54 ± 0.06" 
0.55 ± 0.04 



S Eridani 



5050 = 
2.2 ± 
110.58 
9.0 4 
3.52 ± 

3.75: 



: 100" 

0.9'' 
tO.88" 
0.1 
0.02" 
tO.2 



' iDe Ridderetai]j2006tJ) . '' lTavlodil999l) . " ICarrier et aljj2003h . '' iGlebocki & StawikowskiH200d) . " lESAHl997h 
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Fig. 2. Pow er spectrum of 77 Ser pentis. The lower black one is 
obtained bv^ Carrier et alJ ( l2006l) . the grey one in the middle is 
obtained from (v) and the top one is obtained from the bisector 
velocity span. For clarity the latter two are shifted. 



cross correlation profile, which has an increased signal to noise 
ratio: 



SNR, 



cross-cor ' 



(1) 



with M the number of spectral lines used for the cross correla- 
tion and < iSNKf > the average of the squared signal to noise 
ratios of each spectral line used. The cross correlation profiles 
are constructed from a mask in such a way that the weak, strong 
and blended lines are excluded. They are a good representation 
of an average spectral line of the star whenever the lines are 
formed at not too different line formation regions in the atmo- 

sp here. 

iMathias & Aertsl ( 1 19961) and IChadid et alJ ( 1200 ll) already 
performed a line profile analysis on cro ss correlation pro - 
files for a 6 Scuti star using moments dAerts etalJll992l) . 
Compared to the analysis of a single spectral line this only 
introduces some additional terms in the moments for which 
can be corrected. Furthermore, the line bisect or analysis is 
widely used on cross correlation profiles (see e.g.'Oueloz et al.l 
(2001a), Setiawan et al. (2003), Martinez Fiorenzano et al. 
(120051) . iDall et alJ (I2006t) '). iDall et al.1 (120061) showed that the 
cross correlation bisector contains the same information as sin- 
gle line bisectors. 



Fig. 3. Po wer spectrum o f f Hyd rae. The lower black one is ob- 
tained by Frandsen et alJ (l2002h . the grey one is obtained from 
(v). The two top power spectra are obtained from the bisec- 
tor velocity span. The top one is the original, while the lower 
one is corrected for the led' (11.57yuHz) frequency due to 
the diurnal cycle. For clarity the latter three power spectra are 
shifted. 
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Fig. 4. Pow er spectrum of 6 Er idani. The lower black one is 
obtained bv lCarrieT et al.l(l2003h . the grey one in the middle is 
obtained from (v). The top power spectrum is obtained from the 
bisector velocity span. For clarity the latter two power spectra 
are shifted. 
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The red (sub)giants analysed in the present work show low 
amplitude radial velocity variations, which result in very small 
Une profile variations. The INTER-TACOS (INTERpreter 
for the Treatment, the Analysis and the Correlation of 
Spectra) software pa ckage developed at Geneva Observatory 
(iBaranne et alJl99 a). was used to cross correlate the observed 
spectrum with a mask containing atmospheric lines of red 
(sub)giants with a resolution of 100 ms '. On top of that we 
interpolated to a resolution of 10 ms ' to reach the highest 
possible precision in the computation of the moments. 



2.3. Frequency analysis 

In order to see whether the line moments and line bisector are 
sensitive enough to determine the very small line profile vari- 
ations present in pulsating red (sub)giants, a frequency analy- 
sis is perfo rmed and comp a red wi t h the frequencies obtained 
nrevionslv. IPrandsen et al J (Eioi. lOe Ridderet alJ ^06h), 
ICarrier et al.l ilOOff) and lCarrier et al.l '2003) used radial veloc- 
ity variati ons, obtained with th e optimum weight method de- 
scribed bv lBouchv et alJ(l200ll) . for their frequency analysis. 



2.3.1 . Moments 



A Uri e profile can be described by its moments jAerts et aD 
I1992I) . The first moment (v) represents the centroid velocity 
of the line profile, the second moment (v^) the width of the 
line profile and the third moment (v^) the skewness of the line 
profile. The quantity (v) is thus a particular measure of the 
radial velocity and should therefore show similar frequency 
behaviour. We expect it to perform less well than the radial 
velocity m easure derived from the optimum weight method 
llBouchy et al...2001l) (see below), but (v) is a mode identifica- 
tion diagnostic jAerts et al.lll992l) . while other radial velocity 
measures are not. This is why we re-compute and re-analyse 
the radial velocity from (v) here. 

Frequencies are determined with the conventional method 
of iterativ e sinewaye fit ting ('prewhitening'). A Scargle peri- 
odogram (•Sc argle"l982) is used with frequencies between 
and 15 cycles per day (cd"') (173.6 yuHz) and with a fre- 
quency step of 0.0001 cd ' (0.001 yuHz) for the red giants. 
For the subgiant frequencies between and 150 cd ' and a 
frequency step of 0.001 cd ' are used. The significance of 
the frequencies is calculated with respect to the average am- 
plitude of the periodogram after prewhitening as described in 
llCuschnig et al. ( 1997). The error s in the frequenc y are obtained 
with the method described by Bregeretal 1 11 19991) " and have typ- 
ical values of 10"^ cd"' (lO^V^z) for the giants and 10"^ c d"' 
(lO^'^yuHz) for the subgiant. 

The frequencies of (v) are compared to the frequencies ob- 
tained from radial velocities derived from the optimum weight 
method, mentioned in earlier publications. In Figs m2||3| and 
l^the power spectra of e Ophiuchi, 77 Serpentis, ^ Hydrae and 6 
Eridani are shown, resp ectively. The lower black power spectra 
are the ones obtained by'De Ridder et al. 'fe006bl).ICarrie ret alJ 
(2006) . Frandsen et al. (2002) and Carrier et al. (2003) . while 



the grey ones in the middle are obtained from (v) in the present 
work. 

Although the same data is used, a comparison is made be- 
tween power spectra of two diff'erent representations of the ra- 
dial velocity variations, computed with different methods. The 
radial velocity derived from the optimum weight method is in- 
trinsically more precise than the one from (v) of the cross corre- 
lation. It not only uses the full spectrum instead of a box shaped 
mask, but also relies on weights of the individual spectra. These 
diff'erences influence the noise and the amplitudes in the power 
spectra. On top of that the reference point of each night is deter- 
mined by the observation with the highest signal to noise ratio 
in case of the optimum weight method and by the average of 
the night in case (v) is used. As nightly reference points are 
eff'ectively a high pass filter, different filters induce differences 
at low frequencies in the power spectra. 

Despite these differences the comparison between the 
power spectra in the literature and obtained from (v) is very 
good, especially for e Ophiuchi and 6 Eridan i. For rj Serpentis, 
the dominant frequencies obtained bv lCarri er et al. (2006) are 
slightly different. This is mainly due to the additional ELODIE 
data, which contains a substantial number of observations and 
therefore alters the time sampling of the data. In case only the 
CORALIE data are taken into account, the dominant frequen- 
cies match very well. For ^ Hydrae the dominant frequencies 
obtained by Frandsen et al. (2002) are more distinct than those 
obtained from (v). This is due to the low amplitude of this star, 
which implies a larger relative dilTerence between the differ- 
ent computations of the radial velocity than for the other two 
giants. Nevertheless, the overall shape of the power spectra ob- 
tained from the radial velocity and from (v) is comparable. The 
fact that (v) shows the same behaviour as the radial velocity 
obtained with the method described by iBouchv et alJ ( 1200 ll) . 
indicates that this moment diagnostic is able to detect low am- 
plitude oscillations of red (sub)giants. 

The behaviour of (v^) obtained from a cross correlation is 
not ex actly the same as for a single spectral line. Ch adid et alJ 
EM?) show that the constant of a sinusoidal fit through (v^) 
is different in case of a cross correlation compared to a single 
spectral line. Furthermore, a fit through (y-*) obtained from a 
cross correlation contains a non zero constant, which is not the 
case for a single spectral line. This indicates that we should 
be very cautious in interpreting the constants of the moments. 
Chadid et al. (2001) attribute the different behaviour of (v^) 
and (v^) obtained from cross correlation to the influence of 
possible blending with very faint lines which alters the abso- 
lute values of (v^) and (y-'), but not their variation. 

For the observed red giants, (v^) shows variations in the 
average value per night, which is probably caused by chang- 
ing instrumental conditions. Although a correction for this be- 
haviour is applied by shifting the values of each night to the 
average value of all observations, (v^) does not behave as pre- 
dicted from theory. Frequencies for (v^) are expected to be: 
V,-, Vj + Vj or Vj - Vj , with v,j- the frequencies obtained in (v) 



jMathias et alJl994^ . but are not found in the observations. For 
all four stars, the frequencies obtained for (y-') are the same as 
for (v), which is as expected. 
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2.3.2. Line Bisector 



3. Theoretical mode diagnostics 



The line bisector is a measure of the displacement of the centre 
of the red and blue wing from the core of the spectral line at 
each residual flux. The line bisector is often characterised by 
the bisector velocity span, which is defined as the horizontal 
distance between the bisector positions at fractional flux levels 
in the top and botto m part of the spectral line, see for instance 
iBrown et al.l (Il998 h. The bisector velocity span is supposed to 
remain constant over time in case of substellar companions. 

The bisector velocity span is calculated in this work as the 
difference between the bisector at a fractional flux level of 0.80 
and 0.20. The frequencies of the bisector velocity span are de- 
termined in the same way as the frequencies of the moments. 

Bisector velocity spans are calculated for simulated spectral 
lines of stars pulsating with different modes, using e Ophiuchi's 
amplitude. In case of noiseless line profiles and infinite mode 
lifetimes, the dominant frequency of modes with m = is 2v, 
with V the dominant frequency obtained for (v). In case m 0, 
V is obtained as the dominant frequency. For infinite mode life- 
times the dominant frequency obtained for modes with m = is 
low, but at 2v the power spectrum also shows a clear excess. For 
modes with m it we obtained v as the dominant frequency. 
This is completel y in line with the behaviour expected for (v^) 
kerts et al]l992l) . 

To give an estimate of the minimum signal to noise ratio 
needed to detect a dominant frequency, with the amplitude de- 
tected for e Ophiuchi in the bisector velocity span, we added 
white noise to the simulated spectral lines. These simulations 
reveal that for pulsations with infinite as well as finite mode 
lifetimes with m = 0, a signal to noise ratio of 100000 is 
not enough to reveal a frequency in the bisector velocity span, 
while a signal to noise ratio of order a few 10 000 would suffice 
for modes with m + 0. The signal to noise ratio of the cross 
correlation profiles of our targets does not exceed a few 1 000 
(see Equation [0. This signal to noise level is much less than 
the minimum value needed to detect the oscillations. 

The top graphs in Figs HE] and 0] show the power spec- 
tra of the bisector velocity span for e Ophiuchi, r\ Serpentis, ^ 
Hydrae and 6 Eridani, respectively. The bisector velocity spans 
for e Ophiuchi and r] Serpentis do not show any dominant fre- 
quency. The power spectrum of the bisector velocity span of 
^ Hydrae shows peaks at 1 cd ' (ll.STyuHz) and at its 1 day 
aliases (Fig.|3ltop). The aliases are due to the diurnal cycle and 
in case these are removed no dominant frequencies in the bisec- 
tor velocity span can be seen (Fig. |3 second power spectrum 
from the top). For d Eridani the bisector velocity span shows 
some peaks at low frequencies, but from asteroseismological 
considerations solar like pulsations can not occur at these low 
frequencies in subgiants. These peaks are probably caused by 
instrumental effects. From the present test case, one can thus 
conclude that, at least for low amplitude oscillations, discrim- 
ination between exoplanet companions and pulsations is not 
possible with the bisector velocity span on data wi th reahstic 
signal to noise ratio. This result is consistent with Pall et all 
ll2006l) . 



In order to characterise wavenumbers (l,m) for the modes 
present in the red ( sub)giants, the mat hematical description 
of the line bisector jBrown et al.' ll998h and of the moments 
dAerts et alJl99 2|) are i nvestigated. 

iBrown et al.l (Il998h introduced a mathematical description 
for the line bisector in terms of orthogonal Hermite func- 
tions. The Hermite coefficients hi describe the line shapes. 
Information on the oscillation mode, incl ination and velocity 
parameters of the star was obtained by 'Brown et alJ (Il998h 
from a comparison with theoretically generated line profile 
variations during a pulsation cycle, assuming infinite mode 
lifetimes. To obtain the wavenumbers of the oscillations, the 
Hermite coefficients or bisector velocity span obtained from 
observations are compared with the ones calculated for the sim- 
ulated line profiles. The main drawback of this description in 
Hermite functions is the lack of a pulsation theory directly cou- 
pled to the Hermite coefficients or bisector velocity span. For 
this reason we do not use it here. 

The variations in the observed moments of the line pro- 
files are compared with their theoretical expectations derived 
from oscillation theory (< Aerts et alJll992h . The moments are a 
function of wavenumbers t, m, inchnation pulsation velocity 
amplitude Up, projected rotational velocity usini and intrinsic 
width of the line profile cr. The comparison is performed ob- 
Jectively with a discrim inant (.Aerts et al. ( 1992) . Aerts ( 1996) . 
briquet & AertsI (120031) ') which selects the most hkely set of pa- 
rameters ({, m, i. Up, V sin /, cr). Due to the fact that several com- 
binations of the wavenumbers and velocity parameters result in 
almost the same line profile variation, the "few" best solutions 
resulting from the discriminant should be investigated carefully 
before drawing conclusions about the mode identity. 

The fact that (v) and (v^*) are sensitive to the very small 
line profile variations and the direct connection between the 
moments and pulsation theory makes moments, in principle, 
suitable for the analysis of pulsating red (sub)giants, provided 
that we test its robustness against small amplitudes and the fi- 
nite lifetimes of the modes. This is what we have done in the 
present work. 

3.1. Discriminant 

In the present analysis the discriminant is determined for £ - 
0, 1,2, inclination angles ranging from 5° until 85° with steps 
of 1°, projected rotational velocity vsini and intrinsic width 
of the line profile cr between and 5 kms"' with steps of 0.5 
km s~ '. Furthermore a limbdarkening coefficient u{A) of 0.5 is 
used jvan Hammell 993). 

Due to the fact that the red (sub)giants show low ampli- 
tude variations, the discriminant might give inclination angles 
close to an inchnati on angle of complete cancellation (lACC) 
dChadid etalJl200lh . Complete cancellation occurs when the 
star is looked upon in such way that the contribution to (v) of 
each point on the stellar disk exactly cancels out the contribu- 
tion of another point on the stellar disk. Therefore the angles 
close to 0° and 90° are not taken into account, as well as solu- 
tions close to the lACC at / - 54.7° for i - 2, m - 0. 
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As the discriminant is designed for oscillations inducing 
moment variations larger than the equivalent width of the hne 
(moment of order zero), a check needs to be performed whether 
this technique is applicable in the present work. Therefore three 
series of spectral line profiles with infinite lifetimes and differ- 
ent pulsation velocity amplitudes are generated using the 628 
observation times of e Ophiuchi. The dominant frequency of e 
Ophiuchi, 5.03 c d ' (58.2 //Hz), is used as an input parameter 
at an inclination of / - 35° for wavenumbers £ - 0, 1,2 with 
positive m values, a projected rotational velocity ysin/ = 3.5 
kms ' and an intrinsic width cr - A kms ^ Furthermore, the 
equivalent width of the spectral lines is taken to be 10 kms 
while the minimum amplitude of the pulsation velocity {vp) 
is 0.04 kms ' and increases with a factor of 10 for the dif- 
ferent series. For these series the moments are calculated as 
well as the discriminant. For a pulsation velocity amplitude 
of 0.04 kms which resembles the data best, (v^) of modes 
with m - all show the expected harmonic of the frequency 
(Mathi as et al.lll994l) . but with a different constant. The dis- 
criminant values only differ by an order of 0.01 km s ' . This is 
also the case for a pulsation velocity amplitude of 0.4 kms"', 
which indicates that the discriminant is not suitable for the 
analysis of small line profile variations with an amplitude be- 
low 1 kms"'. For a pulsation velocity amplitude of 4.0 kms"' 
the difference between the discriminant values increases and 
the correct mode is present among the best options. This indi- 
cates that the discriminant is a useful analysis tool for spectral 
line variations with amplitude larger than about half the equiv- 
alent width of the spectral line, but not for red (sub)giants. 



3.2. Amplitude and phase distribution 

As the discriminant is not useful in the present analysis, the 
spectral lines are investigated by comparing the amplitude and 
phase across the line profile with the ones obtain ed from sim- 
ulations. This is the same p rocedure as used by iTelting et alJ 
( Il997h for /3 Cephei and bv lOe Cat et all 12005) for the com- 
plex line profile behaviour due to the multi periodic gravity 
mode oscillations in slowly pulsating B stars. The amplitude 
and phase for each velocity in the cross correlation profile are 
determined by fitting a harmonic function, with the dominant 
frequencies of (v), to the flux values at each velo city pixel of 
the time series of spectra I Schriivers et aljfl997l) . The ampli- 
tudes and phases as function of velocity determined from sim- 
ulated line profiles with £ = 0, 1 , 2, and positive m values are 
shown in Fig. |5] assuming infinite lifetimes. The simulations 
are shown for / = 35°, but the shape of the amplitude and phase 
distribution does not change significantly with inclination an- 
gle. 

The amplitude and phase distribution across the line pro- 
file of € Ophiuchi, 77 Serpentis, ^ Hydrae and S Eridani are 
shown in Figs |6l |H1 and |5] respectively. The dominant fre- 
quencies obtained from (v) are used for the harmonic fits at 
each velocity in the line profile. The frequencies obtained for 
each star are listed in Tables|201E]and|5] The windowfunction 
for 7] Serpentis changes slightly in case only CORALIE data is 
used compared to the one obtained from both CORALIE and 



Table 2. Comparison b etween frequenc i es fouii d in (v) (v(v>) 
and the ones found bv IPe Ridder et al ] ll2006bl) in the radial 
velocity of e Ophiuchi. The significance of V(v> is calculated 
with respect to the average amp litude of the periodogr am af- 
ter prewhitening as described bv lKuschnig et alJ ( ll997l) : 4cr 
99.9% confidence interval, 3.6cr oc 99% confidence interval 
and 3. 25 cr oc 95% confidence interval. 







significance 


v° 


v" 


cd"' 


fiHz 




cd"' 




5.03 


58.2 


4.690- 


5.03 


58.2 


5.46 


63.2 


3.990- 


5.44 


63.0 


5.83 


67.5 


3.980- 


5.82 


67.4 








4.59 


53.1 


3.17 


36.7 


3.430- 


4.17 


48.3 








4.46 


51.6 








5.59 


64.7 








6.19 


71.7 








5.18 


59.9 








6.00 


69.5 








6.43 


74.4 



iDe Ridder et"ai]j2006lJ) 



Table 3. Comparis on between f requencies found in (v> and 
the ones found by Carrier et alJ (I2OO6 ) in the radial velocity 
of 7] Serpentis. The significance (sign) of v<v) is calculated 
with respect to the average amplitude of the per iodogram af- 
ter prewhitening as described by K uschnig et al.l (il997.:) : Act oc 
99.9% confidence interval, 3.6cr oc 99% confidence interval 
and 3.25cr oc 95% confidence interval. 



^<v> 


>'<v> 


signi- 


y" 


y" 




y^ 


cd"' 




ficance 


cd"' 


yUHZ 


cd"' 




11.71 


135.5 


3.8O0- 


10.74 


124.3 


11.74 


135.9 


10.38 


120.2 


3.6I0- 


10.38 


120.2 


10.40 


120.4 












11.49 


133.0 












13.27 


153.6 












10.96 


126.8 












10.90 


126.2 


11.17 


129.3 


4.230- 


11.17 


129.3 


10.33 


119.6 












11.69 


135.3 












7.48 


86.6 












7.28 


84.3 












6.20 


71.8 



" Only CORALIE data from ICarrier et alj tOOdl) taken into 

count; 

* ICarrieretalJj2006h 



ELODIE data jCarrier et al.l2006l) . For the CORALIE data we 
find an extra aliasfrequency of 0. 17 c d"' . The frequency 11.17 
c d"' , only obtained for the CORALIE data, is therefore an alias 
of the 10.33 cd"' frequency obtained by Carrier et al. (200^. 
The frequency v = 3.00 c d"' for 6 Eridani is indicated in bold 
because this is an alias of the diurnal cycle of the observations 
and not due to solar like oscillations. 
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Fig. 5. The amplitude (left) and phase (right) distributions for simulated Une profiles with infinite damping time at an incUnation 
angle of ; = 35°. Top left: ^ = 0, m = 0, top right: { = I, m = 0, middle left: ^ = 1, m = 1, middle right: { = 2,m = 0, bottom 
left: { -2,m- 1, bottom right: € = 2,m = 2 




velocity [km/s] 




velocity [km/s] 




il, [km/,] 



Fig. 6. The amplitude (top) and phase (bottom) distributions as a function of velocity across the hne profile of e Ophiuchi for 
the three significant frequencies obtained from (v): V(v> = 5.03 c d ' (58.2/zHz) (left) V(v> = 5.46 c d ' (63.2juHz) (middle) and 
v<v> =5.83 c d"^ (67.5yuHz) (right). The mean radial velocity of the star, is found to be approximately -9.4 km s"' . 



4. Simulations 

In order to test the robustness of moment diagnostics and the 
amplitude and phase behaviour across the line profiles against 
finite lifetimes of oscillation modes, spectral lines are gen- 
erated with damped and re-excited modes. The moments of 
these lines are obtained and their dominant frequencies derived. 
Furthermore, the ampUtude and phase distribution across the 
spectral line is investigated. 

The simulations are done for a single lineforming region. 
Therefore, the simulated phase distributions are not necessar- 



ily comparable to the observed ones, which are averaged phase 
distributions over a whole range of line forming regions be- 
cause they are based on a cross correlation function. Our in- 
terpretation will therefore mainly be based on the amplitude 
distribution which represents an average amplitude across the 
whole Uneforming region considered in the cross correlation 
function. 
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Fig. 7. The amplitude (top) and phase (bottom) distributions as a function of velocity across the line profile of t] Serpentis for 
the significant frequencies obtained from (v): v^v) - 11-17 cd"' (129.3juHz) (left), v^v) = 11.71 cd"' (135.5juHz) (middle) and 
v<v) - 10.38 c d"' (120.2fjHz) (right). The mean radial velocity of the star, is found to be approximately 9.4 km s"' . 
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Fig. 9. The amplitude (top) and phase (bottom) distributions as a function of velocity across the line profile of 6 Eridani for 
the significant frequencies obtained from (v): v<v> - 59.40 cd"' (687.5juHz) (left), v^v) = 61.08 cd"' (706.9/vHz) (middle) and 
v<v) = 52.68 c d"' (609.7/iHz) (right). The mean radial velocity of the star, is found to be approximately -6.3 kms"'. 



4.1. Damping and re-excitation equations 

A damped and re-excited oscillation mode is damped by a fac- 
tor e^''', with 77 the damping rate, and re-excited before it is 
able to damp out. As a consequence both the amplitude and the 
phase of the oscillation are time dependent: 

/(f) = A(f) sin(27rvf + ifr{t)). (2) 

To simulate such an o scillator we follow the description of 
IPe Ridder et alJ ll2006ah . and we compute 

fit) = B(t) sin(2;rvf) + C(t) cos(2nvt), (3) 



where we let the amplitudes B and C vary with a first order 
autoregressive process in a discrete time domain with time step 
Af: 

B„ = e-''^'B„-i + e„+i (4) 

where £„+i is a Gaussian distributed excitation kick. For more 
details we refer to De Ridder et al. (2006a). The above is ap- 
plied to the equations for the pulsation velocity which are ob- 
tained by taking the time derivative of the displacement com- 
ponents mentioned in Section 3.2, Eq.(2) of De Ridder e t alJ 
(I20n2h . Whenever rotation is neglected, the three spherical 
components of the pulsation velocity, i.e. v,-, vg and v^p, for the 
damped and re-exited case become: 
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Fig. 8. The amplitude (top) and phase (bottom) distributions as 
a function of velocity across the line profile of ^ Hydrae for the 
most dominant frequency obtained from (v): v<v> = 8.42 cd"' 
(97.5yuHz). The mean radial velocity of the star, is found to be 
approximately -4.6 km s"' . 

Table 4. Comp arison between freque ncies found in (v) and the 
ones found by iFrandsen et alJ (l2002l) in the radial velocity of 
^ Hydrae. The significance of v<v> is calculated with respect to 
the average am plitude of the periodo gram after prewhitening 
as described bv lKuschnig et a A significance of 3cr oc 

85% confidence interval. 



Table 5. Comparis on between frequenc ies found in (v) and 
the ones found by ICarrier et alJ ( l2003l) in the radial veloc- 
ity of 6 Eridani. The significance (sign) of V(y) is calculated 
with respect to the average amp litude of the periodogr am af- 
ter prewhitening as described bv lKuschnig et alJ ( ll997l) : Act oc 
99.9% confidence interval, 3.6cr oc 99% confidence interval 
and 3.25cr oc 95% confidence interval. 







significance 
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64.55 


747.1 


65.28 


755.6 


3.7I0- 


66.28 
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" ICarrieretalJj2003h 
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significance 
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y" 


cd-1 
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cd-i 


^Hz 








5.1344(26) 


59.43 








6.8366(27) 


79.13 


8.42 


97.5 
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7.4265(29) 
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8.2318(32) 


95.28 








9.3507(33) 


108.22 
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101.16 








10.0287(43) 


116.07 
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105.13 








8.5339(40) 
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" lFrandseneta'i]<2002h 



s-'1^'-"^'i'-'-':\B„ sm{m(fi + 2nvt) + C„ cos(m^ + 27rvf)), 



Vg 



(5) 



(cos e)) * 

1 

^-n(t-nAt,i,t)^g^^ sin(m^ + Invf) + C„ cos(m<p + Invt)), (6) 



-niKvpN'f 



1 



sinf 



-P'"(cos 6) : 



e 



-il(t-nhtt,ct) 



(Bn cos(nnp + 2nvt) - Cn ?>m{iinp + 2nvt)), (7) 



with Vp proportional to the pulsation amplitude, NJ the nor- 
malisation factor for the spherical harmonics Y"'{0,(p) = 
f""(cos 0)e""'^ and K the ratio of the horizontal to the vertical 
velocity amplitude. 

Line profiles with the same parameters as described in sec- 
tion 3.1 are generated with a d amping time (??'') of two days 
(the estimate for ^ Hydrae by IStelloet alJ (l2004 '). The used 
pulsation amplitude is 0.04 kms The moments, with the 
dominant frequencies and the amplitude and phase distribution 
across the line profiles are investigated. 



4.2. Frequencies 

(v), (v-) and (v^) are determined for the generated series of 
line profiles and the frequencies are obtained in the same way 
as described in section 2.3.1. Rather than fitting a Lorenz pro- 
file to the power spectrum, we simply used the highest peak 
to obtain the oscillation frequency. This implies that the ob- 
tained frequencies diff'er slightly from the input frequency as 
predicted by the Lorenzian probability distribution. The fre- 
quencies obtained from (v-') are the same as the ones obtained 
for (v), which is as expected from theory and also seen in the 
observations of the four red (sub)giants. None of the modes re- 
veal a relation between the frequencies obtained from (v) and 
(v^), which is due to the small amplitudes as described in sec- 
tion 2.3.1. 
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Fig. 10. The amplitude (left) and phase (right) distributions for simulated hne profiles with a two day damping time at an incli- 
nation angle of / - 35°. For each mode ten different realisations are shown in grey, with the average shown in black. Top left: 
£ - 0,m - 0, top right: { - l,m - 0, middle left: £ - l,m - 1, middle right: f - 2, m - 0, bottom left: f - 2, m - 1, bottom 
right: { - 2, m - 2 



4.3. Amplitude and phase distribution 

The amplitude and phase distribution as a function of the ve- 
locity across the line profiles are investigated for the simulated 
spectral lines. The profiles for the cases with ( = 0, 1,2 and 
positive m values, with a damping time of ?/ = 2 days, are 
shown in Fig. ^| The dominant frequency obtained from (v) 
is used for the harmonic fits. This frequency is different for 
each realisation due to the stochastic nature of the forcing. For 
each mode ten realisations are shown in grey, with the aver- 
age shown in black. The shape of the amplitude of axisymmet- 
ric and tesseral modes do not change, although the amplitudes 
per velocity pixel can have different values for different reali- 
sations. The drop to zero amplitude in the centre of the line is 
very characteristic for all the amplitude distributions. However, 
for the sectoral modes, the shape and value of the amplitudes 
may be different for different realisations and the amplitude dis- 
tribution does not necessarily drop in the centre of the line for 
such modes. The behaviour does not change significantly with 
inclination. 

5. Interpretation 

The amplitude distributi on across the line p rofil e is used for 
mode identification, as in lTeltin^ et alJ ll 19971) and lPe Cat et al.l 
For each star we find different distributions for different 
dominant mode frequencies. This clearly indicates that modes 
with different wavenumbers {{, m) must be present in the data. 
Thus, we detect non-radial modes in the observed cross corre- 
lation profiles. In order to identify the wavenumbers of the indi- 
vidual modes, observed amplitude distributions are compared 
with simulated ones, as is commonly done for mode identifica- 



tion. Only very recently, 'Zima" '("2006^ implemented a statistical 
measure to quantify the identification of {{, m) as a function of 
amplitude across the profile for modes with infinite lifetime. 
A similar quantitative measure in the case of damped oscilla- 
tions is not yet available, but visual inspection of amplitude 
diagrams allows a clear discrimination between axisymmetric 
and non-axisymmetric modes which is what we do here, based 
on simulations. The averaged amplitudes over different realisa- 
tions (indicated in black in Fig. I10> are used for this compar- 
ison, keeping in mind that the observations are single realisa- 
tions. Due to the fact that single realisations of different modes 
are in some cases comparable, different sets of wavenumbers 
are occasionally likely for observed modes. 

For € Ophiuchi the amplitude of the first frequency v = 5.03 
cd"' (58.2//Hz) resembles the one of modes with { - 2, m = 
2 mode. The mode with frequency v = 5.46 cd ' (63.2juHz) 
resembles an m - 1 mode, although it does not exclude an m = 
mode. The mode with v - 5.83 c d ' (67.5/iHz) resembles an 
( - 2, m - 2 mode. 

The modes with frequency v - 11.17 c d ' (129.3/zHz) and 
V - 1 1 .71 c d"' (135.5yuHz) obtained for rj Serpentis resembles 
a realisation of the { - 2, m - 2 mode, while the mode with the 
third frequency (v = 10.38 c d ' (102.2/iHz)) is likely to have 
m - Q. 

For ^ Hydrae the most dominant frequency v = 8.42 cd ' 
(97.5//Hz) reveals a mode with either m = or m = 1. 

For 6 Eridani the dominant mode (v - 59.40 cd"' 
(687.5yuHz)) resembles an m - 1, while the second one (v - 
61.08 cd"' (706.9juHz)) is characteristic for an m - mode. 
The third mode (v = 52.68 cd"' (609.7juHz)) shows an addi- 
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tional dip in the second peak, but w = 1 as well as m = seem 
to be possible. 

From the comparison of the observations with the simula- 
tions we can conclude that we see modes with m + and thus 
non-radial modes among the dominant modes in the (sub) gi- 
ants. This result is robust against a change in the value of the 
inclination angle for the simulations. 
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6. Discussion and conclusions 

We have implemented a line profile generation code for 
damped and re-excited solar like oscillations. We found that the 
quantities (v), (v-') and amplitude across the profile are good 
line diagnostics to characterise these oscillations. Our simula- 
tions were made for a single lineforming region while we com- 
pared them with cross correlation profiles based on a mask. For 
this reason, the phase variations across the profile are difficult 
to interpret, since they are an average over the depth in the at- 
mosphere, and this average depends on the shape of the eigen- 
functions. Moreover, we learned from our simulations that the 
phase behaviour can be quite different for realisations, as illus- 
trated in Fig.^l More extensive simulations across the depth 
would be needed to interpret the observed phase behaviour in 
detail. In order to do so, we would need to know the shape of 
the excited eigenmodes as a function of depth in the outer re- 
gions of the star 

From the frequency analysis performed on the bisector ve- 
locity span, it becomes clear that it is not a good diagnostic to 
unravel solar Uke oscillations. 

Individual mode identification is not yet possible with the 
method used in the present work. This is not surprising be- 
cause we did a first exploration of how well-working methods 
for modes with infinite lifetimes could be adapted for damped 
oscillations. The discriminant of the moments is not useful at 
the low amplitudes present in the sub(giants) investigated in 
the present work. Nevertheless, from the resemblance between 
the observed amplitude distribution with the simulated ones, 
and comparison between different modes of the same star, we 
presented clear observational evidence for modes with differ- 
ent m-values in the same star for three of our four examples, 
i.e. non-radial modes are present in red (sub)giants. The ampli- 
tudes for ^ Hydrae are too low to make firm conclusions, and 
since we see only one frequency in the line profile variations, a 
comparison of the amplitude shapes for different modes is not 
possible. Therefore we can not exclude that this star has only 
radial modes as proposed by Frandsen et al. (2002). 

iDziembowski et alJ (1200 ll) and iHoudek & Goughl (l2002h 
show theoretically that stochastically excited radial pulsation 
modes are most likely to be observed in red giants. They pre- 
dict that, for non-radial modes, damping effects will have in- 
fluence in the p-mode cavity as well as in the g-mode cavity, 
while for radial modes only damping in the p-mode cavity is 
present. Therefore, according to this theory, radial modes are 
less severely influenced by damping than non-radial modes. 
Our result pointing towards the presence of non-radial modes, 
calls for a re-evaluation of the theoretical predictions. 
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